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We present evidence for the 6 —* d penguin-dominated decays B + — + K°K + and B° — > K°K° 
with significances of 3.5 and 4.5 standard deviations, respectively. The results are based on a sample 
of 227 million T(4S) -> BB decays collected with the BABAR detector at the PEP-II asymmetric- 
energy e + e~ collider at SLAC. We measure the branching fractions B(B + — > K°K + ) — (1.5 ± 
0.5 ± 0.1) x 10~ 6 (< 2.4 x 10~ 6 ) and B(B° -> K°K°) = (l.lQ±%f 5 ± 0.13) x 10~ 6 , where the 
uncertainties are statistical and systematic, respectively, and the upper limit on the branching 
fraction for K°K + is at the 90% confidence level. We also present improved measurements of the 
charge-averag ed branching fraction B(B + -> K°tt + ) = (26.0 ± 1.3 ± 1.0) x 10~ 6 and CP-violating 
charge asymmetry Acp (K°n + ) = —0.09 ± 0.05 ± 0.01, where the uncertainties are statistical and 
systematic, respectively. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



Flavor-changing neutral currents are forbidden at first 
order in the standard model, but can proceed through 
weak interactions that are described by one-loop "pen- 
guin" diagrams. Such decay processes were first estab- 
lished in the B system more than a decade ago through 
observation of the radiative decay B — > K*j [lj, which is 
dominated by the b — > sj electromagnetic-penguin am- 
plitude. Recently, the analogous gluonic-penguin process 
b — > sg {g — > ss) has been used extensively to test the 
standard model predictions for the CP-violating asymme- 
try amplitudes of decay modes such as B° — > 4>K° [2j ■ To 
date, no direct evidence has been found for decays dom- 
inated by the corresponding b — > dg transition, whose 
amplitude is suppressed relative to that for the b — > sg 
process by the small ratio Kd/^ts involving elements 
of the Cabibbo-Kobayashi-Maskawa quark-mixing ma- 
trix 0]. Such decays could play an important comple- 
mentary role in the search for new physics in the B sys- 
tem. 

In this Letter, we report evidence for the decays B + — > 
K°K + and B° — > K°K°, which are expected to be dom- 
inated by the b — > dg (g — > ss) penguin diagram, and 
an updated measurement of the branching fraction and 
direct CP-violating charge asymmetry for B + — > K°tt + 
(the use of charge conjugate modes is implied throughout 
this paper unless otherwise stated). Our previous search 
for the KK° modes yielded branching-fraction upper lim- 
its at the level of 2 x 10 -6 0, which are consistent with 
recent theoretical estimates based on perturbative calcu- 
lations as well as the lower bounds implied by SU(3) 
symmetry 0. 

Once the decay B° — > K°K° has been established, a 
measurement of its time-dependent CP-violating asym- 
metry (through the technique described in Ref. |2j) could 
provide important constraints on physics beyond the 
standard model. Assuming top-quark dominance in the 
penguin loop, the asymmetry is expected to vanish in 
the standard model while contributions from super- 
symmetric particles could be significant 0. Although 



soft rescattering effects could weaken the sensitivity to 
new physics in this mode 0, the ratio of decay rates 
for B + — > K°K + , K°ir + can be used to constrain the 
relative size of such effects fi"l| . 

Recent measurements of the partial-rate asymmetry in 
B° -> K+ir- decays by the BABAR G3 and Belle ex- 
periments have established direct CP violation in the B 
system. In this Letter, we search for direct CP violation 
in the decays B + — » K°tt + , K°K + through measure- 
ment of the charge asymmetry 

. _ r(£r -> r)-r(B+ -> /+) 
r(s- ->/-) + r(s+ -►/+)' 

where /± = X°tt ± , K^K*. The decay B~ -> K°tt- 
is dominated by the b —> s penguin process and, ne- 
glecting rescattering effects [llj, is expected to yield 
Acp ~ 1% HE1- 

Observation of a significant charge 
asymmetry could therefore indicate new physics entering 
the penguin loop |l5j . The decay rate and charge asym- 
metry in K®ir + can also be used to constrain the angle 
7 of the unitarity triangle [l6| . 

The data sample used in this analysis contains (226. 6± 
2.5) x 10 6 T(45) -► BB decays collected by the BABAR 
detector 

.17] 

at the SLAC PEP-II asymmetric-energy 
e + e~ collider. The primary detector elements used in this 
analysis are a charged-particle tracking system consisting 
of a five-layer silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) surrounded by a 1.5-T solenoidal 
magnet, and a dedicated particle-identification system 
consisting of a detector of internally reflected Cherenkov 
light (DIRC). 

We identify two separate event samples corresponding 
to the decay topologies B° ->■ K°K° and B+ -> K°h+ , 
where is either a pion or a kaon. Neutral kaons are re- 
constructed in the mode K° — > tt + tt~ by combining pairs 
of oppositely charged tracks originating from a common 
decay point and having a t7 + tt~ invariant mass within 
11.2 MeV/c 2 of the nominal mass 0|- To reduce com- 
binatorial background, we require the measured proper 
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decay time of the K® to be greater than five times its 
uncertainty. Candidate h + tracks are assigned the pion 
mass and are required to originate from the interaction 
region and to have an associated Cherenkov angle (6 C ) 
measurement with at least six signal photons detected 
in the DIRC. To reduce backgrounds from protons and 
leptons, we require 9 C to be within 4 standard deviations 
(a) of the expectation for either the pion or kaon parti- 
cle hypothesis. The B° sample is formed by combining 
pairs of K® candidates, while the B + sample is formed 
by combining and h + candidates. 

We exploit the unique kinematic and topological fea- 
tures of charmless two-body B decays to suppress the 
dominant background arising from the process e + e~ — > 
QQ (<1 = u,d,8,c). For each B a candidate, we require 
the difference AE between its reconstructed center-of- 
mass (CM) energy and the beam energy (y / s/2) to be 
less than 100 MeV. For B + candidates, we require 
— 115 < AE < 75 MeV, where the lower limit ac- 
counts for an average shift in AE of —45 MeV in the 
K°K + mode due to the assignment of the pion mass 
to the K + . We also define a beam-energy substituted 
mass Toes = \/ (s/2 + p; ■ pB) 2 /Ef — p 2 B , where the B- 
candidate momentum p b and the four- momentum of the 
initial e + e~ state (Ei,pi) are calculated in the labora- 
tory frame. We require 5.20 < toes < 5.29GeV/c 2 for 
B candidates in both samples. To suppress the jet-like 
qq background, we calculate the CM angle 8g between 
the sphericity axis of the B candidate and the sphericity 
axis of the remaining charged and neutral particles in the 
event, and require |cos(0g)| < 0.8. 

After applying all of the above requirements, we find 
1939 (20441) candidates in the B° (B + ) samples, respec- 
tively. The fraction of events containing more than one 
B candidate is negligible (< 0.5%). The total detection 
efficiencies are given in Tableland include the branching 

8| and a probability of 50% 
We use data and simulated 
Monte Carlo samples |2fJ to verify that backgrounds from 
other B decays are negligible. The selected samples are 
therefore assumed to be composed of signal B decays 
and background candidates arising from random combi- 
nations of tracks and mesons in qq events. 

To determine signal yields in each sample, we apply 
separate unbinned maximum-likelihood fits incorporat- 
ing discriminating variables that account for differences 
between BB and qq events. In addition to the kinematic 
variables toes and AE, we include a Fisher discriminant 
T defined as an optimized linear combination of the 
event-shape variables J2iPi and J2iPi cos2 (^*)j where p* 
is the CM momentum of particle i, 9* is the CM angle 
between the momentum of particle i and the i3-candidate 
thrust axis, and the sum is over all particles in the event 
excluding the B daughters. 



The likelihood function to be maximized is defined as 



jv 



fraction for K® — *■ ir + ir 
K s K s JZ ! 



for K°K° 



£ = cx P y-J2 n >J n 



where and Vi are the yield and probability density 
function (PDF) for each component i in the fit, and N 
is the total number of events in the sample. For the B° 
sample there are only two components (signal and back- 
ground), and the total PDF is calculated as the product 
of the individual PDFs for toes, AE, and T . We com- 
bine B + and B~ candidates in a single fit and include 
the PDF for C to determine separate yields and charge 
asymmetries for the two signal components, K®ir and 
KgK, and two corresponding background components. 
For both signal and background, the K^h^ yields are 
parameterized as n± = n(l =F Acp)/2', we fit directly for 
the total yield n and the charge asymmetry Acp- 

The parameterizations of the PDFs are determined 
from data wherever possible. For the B + sample, the 
large signal K®ir + component allows for an accurate de- 
termination of the peak positions for toes and AE, as 
well as the parameters describing the shape of the PDF 
for T . We therefore allow these parameters to vary freely 
in the fit. The remaining shape parameters describing 
toes and AE are determined from simulated Monte Carlo 
samples and are fixed in the fit. Except for the mean 
value of AE, which is shifted by our use of the pion 
mass hypothesis for the h + candidate, we use the K®ir + 
parameters to describe signal K®K + decays. The pa- 
rameters describing the background PDFs in toes and T 
are allowed to vary freely in the fit, while the AE pa- 
rameters are determined in the signal-free region of toes 
(5.20 < toes < 5.26 GeV/c 2 ) and fixed in the fit. For 
both signal and background, the 9 C PDFs are obtained 
from a sample of D* + — ► D°ir + (D° — > K~ir + ) decays 
reconstructed in data, as described in Ref. [13. For the 
B° sample, all shape parameters describing the signal 
PDFs are fixed to the values determined from Monte 
Carlo simulation, while the peak positions for Toes and 
AE are derived from the results of the fit to the B + 
sample. We allow the background T shape parameters 
to vary freely, while the PDF parameters for toes and 
AE are fixed to the values determined from data in the 
signal-free regions 100 < \AE\ < 300 MeV (for to E s) and 
5.20 < toes < 5.26GeV/c 2 (for AE). 

Several cross-checks were performed to validate the fit- 
ting technique before data in the signal region were ex- 
amined. We confirmed the internal self-consistency of the 
fitting algorithm by generating and fitting a large set of 
pseudo-experiments where signal and background events 
were generated randomly from the PDFs with yields cor- 
responding to the expected values based on our previ- 
ous analysis of these modes Q. The fitted signal yields 
for all three modes were unbiased. Correlations among 
the discriminating variables in background data events 
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TABLE I: Summary of results for the total detection efficiencies e, fitted signal yields n, signal-yield significances s (includ- 
ing systematic uncertainty), charge-averaged branching fractions B, and charge asymmetries Acp (including 90% confidence 
intervals). The efficiencies include the branching fraction for — > -k + tx~ - and the probability of 50% for K°K° K^K%. 
Branching fractions are calculated assuming equal rates for T(45) — » B°B° and B + B~ [2^| . For K°K + , we give both the 
central value of the branching fraction and, in parentheses, the 90% confidence-level (CL) upper limit. 



Mode 




e (%) 


n 


»(«7) 


S (HT 6 ) 


Acp 


Acp (90%CL) 


B+ -> 


K°7T+ 


12.6 ± 0.3 


?44 +37 +21 
'**-36 -17 




26.0 ± 1.3 ± 1.0 


-0.09 ±0.05 ±0.01 


[-0.16, -0.02] 


B+ -> 


K°K+ 


12.5 ± 0.3 


a 1+15 +3 
* A -13 -2 


3.5 


1.5±0.5±0.1(< 2.4) 


0.15 ±0.33 ±0.03 


[-0.43,0.68] 


B° — > 


K°K° 


8.5 ±0.6 


23t 8 7 ± 2 


4.5 


1.19i°;«±0.13 







are found to be negligible. To check for residual cor- 
relations between the discriminating variables in signal 
events, we performed a second test for the K®K® mode 
where simulated Monte Carlo samples of signal events 
were mixed with background events generated directly 
from the PDFs. We observed an average bias correspond- 
ing to approximately one event and include this effect in 
the systematic uncertainty on the fitted K®K® yield. Po- 
tential K®ir — > K®K cross-feed was evaluated by fitting 
large samples of simulated Monte Carlo signal events. 
The resulting small (< 0.5%) biases are included in the 
systematic uncertainty on the fitted yields. 

The fit results supersede our previous measurements of 
these quantities and are summarized in Table [I] The sig- 
nal yields for B + — > K° S K + and B° — > K®K® correspond 
to significances of 3.5a and 4.5cr (including systematic 
uncertainties 0), respectively, and are consistent with 
our previous results [j], as well as with the results of 
other experiments j24|. The signal yield for B + — ► AT°7r + 
is somewhat higher than expected from our previous re- 
sult. A re- analysis of the first 88 million BB events yields 
285 ±21 KgTr + signal events, compared with 255 ±20 re- 
ported in Ref. Q . Approximately half of this difference is 
due to reprocessing of the data with improved calibration 
constants. The remaining difference is due to improved 
knowledge of the PDF parameters, which were the largest 
source of systematic uncertainty for the previous result. 
We find 459 ±29 events in the remaining 139 million BB 
events, which is consistent with the signal yield obtained 
in the first part of the sample. 

For the K®K + mode, we compute an upper limit 
on the signal yield as the value of no f° r which 
Jo™ £ max d7j/ J °° £maxdn = 0.9, where £ max is the likeli- 
hood as a function of n, maximized with respect to the re- 
maining free parameters. The corresponding branching- 
fraction upper limit is calculated by increasing n and 
reducing the efficiency by their respective systematic un- 
certainties. 

We compare data and PDFs in the high-statistics 
K®ir + mode using the event-weighting technique de- 
scribed in Ref. 25]. For the plots in Figs.n( a 7 D )j we 
perform a fit excluding the variable being shown; the 



covariance matrix and remaining PDFs are used to de- 
termine a weight that each event is either signal (main 
plot) or background (inset). The resulting distributions 
(points with errors) are normalized to the appropriate 
yield and can be directly compared with the PDFs (solid 
curves) used in the fits. We find good agreement be- 
tween data and the assumed PDF shapes for toes and 
AE. In Figs.n(c-f), we show projections of the K®K + 
and K®K® data obtained by selecting on probability ra- 
tios calculated from the signal and background PDFs (ex- 
cept the variable being plotted) . The solid curves in each 
plot show the fit result after correcting for the efficiency 
of this additional selection. 

Systematic uncertainties on the signal yields are due to 
the imperfect knowledge of the PDF shapes. We evaluate 
this uncertainty by varying the PDF parameters that are 
fixed in the fit within their statistical errors, and by sub- 
stituting different functional forms for the PDF shapes. 
For the charged modes, the largest contribution is due 
to the signal parameterizations for wes (-15 events for 
A'°7r+, t{]r events for K° S K+) and AE (+^ 6 events for 
K^n + , events for K®K + ), while for the neutral mode 
it is due to uncertainty in the background toes shape 
(±0.7 events) and the potential fit bias (±1.4 events). 
The systematic uncertainties on efficiency estimates are 
dominated by the selection on cosOs (2.5%) and the un- 
certainty (1.2% per K^) in reconstruction efficien- 
cies evaluated in a large inclusive sample of K® mesons 
reconstructed in data. For the charge-asymmetry mea- 
surement, we use the background asymmetry to set the 
systematic uncertainty We find background asym- 
metries of -0.005±0.010 and -0.002±0.011 for AT°7r and 
KgK events, respectively. Both results are consistent 
with zero bias, and we assign the statistical uncertainty 
(0.01) as the systematic error on Acp(K°n + ). 

In summary, we find evidence for the decays B + — > 
K°K + and B° — > K°K° with branching fractions on the 
order of 10~ 6 and significances of 3.5c and 4.5<r, respec- 
tively, including systematic uncertainties. These results 
represent direct evidence for the b — > dg penguin-decay 
process. The branching fractions are consistent with re- 
cent theoretical estimates implying that soft reseat- 
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FIG. 1: Distributions of (a) jtibs and (b) AE for signal (main 
plot) and background (inset) B + — > K%n + candidates (points 
with error bars) using the weighting technique described in 
the text. Solid curves represent the corresponding PDFs used 
in the fit. In (c-f) we show projections of mes and AE for 
K° S K + (c,d) and K° S K° S (e,f) decays (points with error bars) 
enhanced in signal decays using additional requirements on 
probability ratios. Solid curves represent the PDF projections 
for the sum of signal and background components, while the 
dotted curve shows the contribution from background only. 



tcring effects may not play an important role in these de- 
cays. We also measure the branching fraction B(B + — > 
K° s ir+) = (26.0 ± 1.3 ± 1.0) x 10~ 6 and the CP-violating 
charge asymmetry A C p(K°tt+) = -0.09 ± 0.05 ± 0.01, 
which are both consistent with previous measurements 
by other experiments |2J, |26j . 
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